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Using the PQSCHA semiclassical method, we evaluate thermodynamic quantities of one-dimensional 
Heisenberg ferro- and antiferromagnets. Since the PQSCHA reduces their evaluation to classical-like 
calculations, we take advantage of Fisher's exact solution to get all results in an almost fully analyt- 
ical way. Explicitly considered here are the specific heat, the correlations length and susceptibility. 
Good agreement with Monte Carlo simulations is found for S > 1 antiferromagnets, showing that 
the relevance of the topological terms and of the Haldane gap is significant only for the lowest spin 
values and temperatures. 



Several applications to condensed matter systjems have 
demonstrated the usefulness of the improvecU'El effec- 
tive potential approach^; its generalized version for non 
standard Hamiltonian, the so called pure- quantum self- 
consistent harmonic approximation (PQSCHApfl, has 
also been successfully applied to different spin systems. 

We consider here the one-dimensional isotropic Heisen- 
berg model Hamiltonian, 



(1) 



with the exchange interaction J > restricted to nearest- 
neighbour sites of a simple cubic d-dimensional lattice; 
the sign — refers to the ferromagnet (FM), -I- to the an- 
tiferromagnet (AFM). The thermodynamic quantities of 
this model were successfully calculated for two-ErLl and 
three-dimensionalH magnets, both of which are charac- 
terized by a ground state that can be obtained pertur- 
batively starting from the classical-like minimum-energy 
configuration. 

In one dimension the situation is largely different and 
the ferro- and antiferromagnets, at variance with the 
classical case where both models are mapped onto the 
classical non linear Schrodinger equation, behave in a 
markedly different way. In fact, ferromagnets do not 
present any peculiarity, being their (ordered) ground 
state, as in higher dimension, the quantum counterpart 
of the classical minimum-energy configuration. The rel- 
evant excitations, both linear and non-linear, which con- 
tributes to the thermodynamic properties and destroy 
the long-range order at any finite temperature, are basi- 
cally the same in the quantum and in the classical case. 
The absence of long-range order at any T ^ Oa implies 
that the linear excitations spectrum should be considered 
only up to wavelengths of the order of the correlation 
length ^(T). 

Quantum effects have a much more apparent impact 
on the qualitative behaviour of antiferromagnets. Its 
ground state cannot be obtained perturbatively from the 



Neel configuration, and the relevant quantum excitations 
are completely different from the classical ones. More- 
over, as firstly suggested by HaldaneE3ilil, integer and 
half-integer spin chains display a qualitatively different 
low-temperature behaviour, which arises from a topolog- 
ical term in the path- integral description of spin systems. 
For half-integer spins, the interference term leads to gap- 
less excitations, while the Haldane gap appears for inte- 
ger 5". However, these effects should rapidly disappear 
for increasing spin value, as the interference becomes less 
destructive and the gap vanishes exponentially. 

Hence, in one dimensional magnets should also exist 
regimes where semiclassical methods can be sensibly ap- 
plied, i.e., when the spin length and/or the temperature 
increase and the classical behaviour is approached. In 
such regimes, the PQSCHA is a good tool for calculating 
thermodynamic quantities and should be very competi- 
tive in comparison with other semiclassical ones, as, for 
instance, the theory based on the use of real-SMce coher- 
ent states that has been recently introducedt3. Indeed, 
by the PQSCHA most calculations can be performed an- 
alytically, and with a full quantum inclusion of the linear 
excitations in wave-vector space. 

The final outcome of the application of the PQSCHA 
to Heisenberg models is that the free energy of the quan- 
tum model described by Eq. (P is given by the the free 
energy of the effective classical Heisenberg Hamiltonian, 



Hcff = ±Ijs^ e^{t) J2 sj-sj+s + Njs^g(t) , 



(2) 



whose thermodynamicj-^roperties are exactly known af- 
ter the work by Fisherll^. In Eq. ^ Sj are unit vectors, 
5 = 5+1/2 plays the role of the 'classical' spin length 
and JS^ that of the overall energy scale: we hence de- 
fine t = ksT/JS"^ as the reduced temperature. The 
temperature- and spin-dependent parameters 9{t) and 
Q{t) account for the effects of the pure-quantum fluctua- 
tions. Their explicit expressions are: 



1 



0'fk 



(3) 
(4) 



It is worthwhile pointing out that the first term oiQ(t) rp.- 
stores the quantum free energy of the linear excitationsO. 
The renormalization coefficient 2?(t) reads 



NS V 



coth/fe- 




(AFM) 
(FM) 



(5) 



and represents indeed the pure-quantum nearest-neigh- 
bour transverse spin fluctuations in self-consistent Gaus- 
sian approximation. This essential renormalization co- 
efficient of the PQSCHA for Heisenberg models, as well 
as the connected quantities 9^{t) and G{t), are global 
parameters, i.e., they take into account the quantum ef- 
fects only on average, so that the details of the exci- 
tation spectrum are smeared out. Furthermore, fk = 
{hujk)l{2kBT) = uJk/2St, where 



zn^ (1 -7fc) 
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(6) 



are the dimensionless spin- wave frequencies, whose renor- 
malization factor (t) = i (j)'^ + \/()^ — 4te / is calcu- 
lated by taking into account only the thermal fluctua- 
tions with |A;| > 7r/^(i), so that e = 1 — 1/^ for the 
onc-dimcnsional system. We point out that the contribu- 
tions to the pure-quantum renormalization coefficient (j^) 
are weighted by the Langevin function coth/fc — fk^j so 
that the major role is played by the high-frequency (short 
wavelength) excitations, which are just those that survive 
in spite of the absence of long-rangeiOfidcr. The PQSCHA 
in the low-coupling approximationO'BQ employed to de- 
rive Eq. (^) neglects contributions of order so that V 
must be small compared to one. 

Using the exact results for the classical one-dimen- 
sional Heisenberg model of Ref. we can thus easily 
write an analytical expression for the free energy per spin 
of the quantum spin chain within the PQSCHA: 
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The other macroscopic thermodynamic quantities, e.g. 
internal energy and specific heat, can be easily obtained 
from this equation by (numerical) derivation, taking care 
of the ^-dependence of Q and 0, which prevents us from 
using directly the expressions of Ref. 13| in evaluating 
such quantities. One can also obtain the spin correlation 
function and the susceptibility by means of the formu- 
las reported in Ref. M, where it is shown also how the 
correlation length ^{t) is simply related to its classical 



counterpart only by the change of the temperature scale 
involved in the renormalization of the exchange constant. 



(8) 



This formula is of remarkable simplicity, especially when 
one notiees that ^ci(i) has a very simple analytical ex- 
pressions. 



ln(cotht"^ - t)] ^ 



(9) 




FIG. 1. Pure-quantum renormalization coefficient 5") 
vs. temperature t = T/JS^, for the Heisenberg FM chain. 
From the uppermost curve, the spin values are S — 1/2, S — 1 
and S = 3/2. 




FIG. 2. Specific heat c{t) of the Heisenberg FM chain for 
S — 1/2, S — 1 and S — 3/2, respectively, starting from the 
lower-most curve. 

In principle, the PQSCHA approach should work well 
for the one-dimensional Heisenberg FM, because: i) 
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their ground state is ordered and is an eigenstate of 
'S'tot = X^i ^i'^ absence of long-range order at any 

finite temperature is due to nonlinear excitations, which 
can be well treated in semiclassical (l/S-expansion) ap- 
proximation; in) the pure-quantum fluctuations related 
to the linear excitations are accounted for through the co- 
efficient Therefore, the only limitation we should 
care of for the FM is a possible too high value of V: we 
shall assume that the condition V < 0.5 must be satisfied 
to get reliable results. As shown in Fig. || this occurs at 
any temperature for S > 3/2. As for lower values of S it 
is apparent that for 5* = 1 the approach can be used in a 
wide range of temperatures starting from t 0.25; on the 
other hand, the strongest quantum case, S = 1/2, can- 
not be well described except for very high temperatures 
t>l. 

Unfortunately, to the best of our knowledge, reference 
data on one-dimensional Heisenberg ferromagnets are not 
available for intermediate spipj\«Jues, and we found only 
rather old ones for S = IfMM. However it could be 
appealing to derive some thermodynamic quantities by 
means of the effective Hamiltonian (|^) and Eq. (^. The 
specific heat of the FM spin chain is shown in Fig. ^ 
for spin values S = 1/2, S = 1 and S = 3/2. Comparing 
with Fig. 0we see why the curve for S — 1/2 is truncated 
at i ~ 1: for higher temperature the behaviour pf|Cl(i) is 
in agreement with the available numerical datallii'lia . 



by any semiclassical approach, so that we do not expect 
that peculiar quantum properties at very low temper- 
ature and low values of the spin could be reproduced 
within the PQSCHA. Instead the PQSCHA is expected 
to give good results when the spin becomes larger and 
larger approaching the classical AFM Heisenberg model. 
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FIG. 4. Specific heat c(t) of the Heisenberg AFM chain 
for 5* = 1/2, 1, 3/2, 2, and oo (classical limit), respectively, 
starting from the lower-most curve. 
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FIG. 3. Pure-quantum renormalization coefficient T>{t, S) 
vs. temperature t = T/JS^, for the Heisenberg AFM chain. 
From the uppermost curve the spin values are 5* — 1/2, 1, 
3/2, and 2. 

Turning to antiferromagnets, the temperature be- 
haviour of X'(t) is shown in Fig. Using the same crite- 
rion as above, we can deduce that the PQSCHA should 
work well for any temperatures for S > 1 while its va- 
lidity is confined to t ^ 0.15 for S = 1/2. However, 
we must recall that the afore-mentioned quantum effects 
which strongly modify the nature of the ground state 
and lead to the Haldane gap, cannot be accounted for 




FIG. 5. Comparison of the specific heat c{t) of the Heisen- 
berg AFM chain, plotted vs.j-T/J — tS?, with numerical 
simulation data for S = l/'M, S = M, S = 3/2Ea and 

s = -M. 

The results for the specific heat plotted in Figs. ^ 
and ^ seem to confirm this prediction for the thermo- 
dynamic quantities. The-comparison with the existing 
quantum Monte GarlotJIlJ and transfer-matrix Renor- 
malization Groupli3 data shows that the agreement im- 
prove for larger spin and can be considered very good 
almost at all temperature already for S — 3/2. 
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ture behaviour can not be, as expected, reproduced. 

The conchision is that the thermodynamic quantities 
of the Heisenberg AFM are not strongly affected by the 
consequences of the Haldane ground state at intermediate 
and high temperatures. SmaU differences appear only at 
the lowest values of the spin and from the point of view 
of the PQSCHA it is difficult to ascertain if they are due 
to the "Haldane effects" or to the rather large value ofV. 
On the other hand, it is curious to note that the Haldane 
conjecture was theoretically proven for large S - where 
the role of the topological term decreases and the gap 
is exponentially vanishing, giving negligible quantitative 
contributions to the thermodynamics. 



FIG. 6. Correlation length ^{t) of the AFM Heisenberg 
chain for S = 1/2, 5 = 1, S = 3/2 and S = 2, respec- 
tively, starting from the lower-most curve; the dashed line is 
the classical result. 
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FIG. 7. Comparison of the uniform susceptibility x(i) of 
the AFM Heisenberg chain with Quantum Monte Carlo sim- 
ulation datatJ. Spin values are 5=1/2 (line and upward 
triangles), S = 1 (line and circles), 5 = 3/2 (line and squares) 
and S — 2 (line and downward triangles), respectively, start- 
ing from the lower-most curve; the dashed line is the classical 
result. 

In Figs. ^ and |^ we show the PQSCHA results for 
the correlation length ^(t) and the uniform susceptibility 
x{t) of the AFM chain. The latter quantity has been 
shown to be particularly sensitive to test the peculiar 
quantum effects related to the Haldane gap; neverthe- 
less. Fig. 1^ clearly shows that a semiclassical approach 
like the PQSCHA compares well with numerical data not 
only where the classical regime is already approached, 
but also in the intermediate temperature range, where 
the curves for low spin values already are clearly different 
from the classical behaviour; only the very low tempera- 
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